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ABSTRACT 


The effect of temperature, environment and heat treatment on 
the fatigue crack propagation rates in a fine grained pipeline steel 
corresponding to A. P.L specifications 5L Grade B and/or 5LX-X42 
were investigated. Constant deflection bending fatigue tests at 640 
c. p.m. were conducted on single edge notched specimens which con- 
tained a 1/16 inch diameter key hole notch as a stress riser. A fatigue 
machine was designed and built specifically for this project. The 


fatigue tests were performed at -50, -10, 21 and 71°C in argon and at 


-50 and 21°C in hydrogen at ambient pressure. One set of specimens 
was heat treated to explore the effect of grain size on crack propaga- 
tion rate. The grain size increased from 0.022 mmto0.026 mm. 

The surface plastic zones and reversed plastic zones were 
photographed to determine if any correlation could be made between 
the size or shape of the zones and the fatigue crack growth rates. The 
fatigue fracture surfaces were studied with both the scanning electron 
microscope and the transmission electron microscope. 

In addition to fatigue tests, static tensile tests at strain rates 
of 0. 2 and 0. 002 in/in/min, cyclic tensile tests and impact tests were 
carried out at the four test temperatures. 

The data were evaluated in terms of the crack propagation 


rates (da/dN) as a function of SK, according to da/dN=C SK™, 
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The values of m varied from a low value of 2. 66 at 71°C to a maximum 
of 6. 81 at -10°C after which it decreased to 4. 49 at -50°C. The values 
of C varied inversely to those of m and it was found that a relation- 
ship of the form m - mpg = A ln C holds for a wide variety of materials. 
The experimentally determined relationships were compared to a 
theoretically derived crack growth law due to Tomkins (54) and a close 
correlation is found if cyclic material properties are used. 

Hydrogen at ambient pressure had no significant effect on crack 
initiation life or initial crack growth rates, however after an incuba- 
tion period crack growth rates an order of magnitude faster occurred. 
A critical SK value of 20-22 ksi Jin was sufficient to cause unstable 
crack growth. Reversed plastic zone sizes and cyclic stress-strain 
data show a much better correlation with fatigue properties as these 


properties are affected by strain ageing. 
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1 INTRODUCTION 


With the recent advances in technology low alloy-low carbon 
steels in the 50 to 80 ksi yield strength range are now a reality (1). 
The properties of these new steels are obtained by a combination of 
grain refinement by means of niobium, vanadium, and rare earth 
additions, plus a change in the steel processing practice such as 
rolling to a lower finishing temperature (Low Temperature Thermal 
Mechanical Treatment, ITTMT) and/or rapid cooling through 
special quenches or coiling procedure (1)(2)(3)(4). These steels have 
found uses in such structures as pipelines, refrigerated storage 
vessels, automobile bumpers, engine mounts and frame members (1), 
Failure of such structures or components usually result in serious 
damage. 

Fatigue is one of the most common failure mechanisms, Even 
though fatigue failures have been recognized and studied for more than 
a century, only recently has structural fatigue been recognized as a 
field of major concern, The development of low cost high-strength 
alloys and improved fabrication and joining techniques have led to the 
production of large structures which must sustain repeated stresses in 
service, Whereas smaller components can be designed to prevent 
crack initiation at an increase in cost, large structures must be designed 
economically, thus the prevention of terminal crack propagation must 


be considered where much shorter fatigue lives are involved. 
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Because these advances have taken place in a short period of 
time much technical data is lacking, In a recent symposium (1) 
representatives of automobile manufacturers pointed out that more 
data on formability, corrosion resistance, cyclic yield, and fatigue 
were required. In another recent paper (6), concerned with metals 
and alloys for low temperature use, it was pointed out that very little 
information was available regarding the effect of low temperature on 
fatigue crack propagation, Clark and Trout (7), have studied the effect 
of temperature (-100 to 75°F), and section size on fatigue crack growth 
rates in a forging grade Ni-Mo-V steel and found that as the tempera- 
ture decreased the crack growth rate decreased for a given stress 
intensity. Rolfe and Munse (8), in a study of fatigue crack propagation 
in mild steel found that as the temperature decreased (78 to -40 TF), 
the fatigue crack growth rate decreased, however these steels do not 
resemble the recently developed alloys therefore a close comparison 
can not be made. 

Because little data is available for these steels a study of 
fatigue crack propagation was initiated, Although fatigue may not be 
the cause of the final failure, it is often fatigue which causes sub- 
critical flaws to grow to a size where failure by some other mechanism 
can take over (6), Fatigue crack growth is affected by many variables, 
but test temperature and environment seem to be most important and 


as a result were chosen as the two main variables of this investigation, 
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2.) BAPE RIMENTAL PROCEDURE ANDtPROGRAM 
2.1 Introduction 

The purpose of this investigation was to study the effect of 
varying service temperature on the fatigue crack propagation rates 
in a fine grained pipeline Steel. The temperatures’ selected were 
-50, -10, 21 and 71°C, thus ranging from a low ambient tempera- 
ture to a temperature range where strain ageing effects are signi- 
ficant. To delineate environment effects tests were either performed 
in dried argon or hydrogen at ambient pressure. 

The following data were recorded from the fatigue tests: total 
fatigue life, number of cycles for crack initiation, crack propagation 
rate as a function of the range of stress intensity, fatigue striation 
Spacings, and the size of the surface and reversed plastic zone. 

In order to analyse and correlate the fatigue data the following 
additional tests were performed at the four different temperatures; 
Charpy impact tests, tensile tests at two strain rates (0, 002 and 0. 2 in/ 
in/min), and cyclic tensile tests. 

22 Test Material 

The test material selected meets A.P.L 5L Grade B and/or 
A. P. I, 51X-X42 specifications. The supply of the 3/8" skelp by 
ST ELC O is gratefully acknowledged, The material was tested 
in the as received condition with the exception of one set of test 


specimens which was given a full anneal by heating to 1750 °F for 15 
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minutes and then furnace cooled. The chemical composition and 
mechanical properties are listed in Table l, 
Zao Miecrostnucture 

As in most materials which have been mechanically worked the 
microstructure has an orientation relationship to the direction of 
working. The microstructure in the various orientations with respect 
to the rolling direction is shown in Figure 1. The notation used to 
describe the orientation with respect to the rolling direction is shown 
schematically in Figure 2 and will be used throughout the thesis, 
The microstructure 90° to the rolling direction is a fairly even dis- 
tribution of pearlite and ferrite. Parallel to the rolling direction 
the microstructure is 'banded'", A 'banded" structure consists of 
alternate layers of ferrite and pearlite. The pearlite is separated by 
two to three grains of ferrite, the average size of which is 0. 0223 mm, 

Manganese sulphide inclusions are present as indicated by 
Figure 2A. A,S, T, M, Designation E45 was used to determine the 
inclusion content (9), The microscopic method was used and it was 
found that inclusions of Type A were present; their size and distri- 
bution density placing this steel below the A, S, T, M. numerical 
designation of 1, thus indicating a clean steel, 
2.4 Charpy Impact Test Results 

Impact tests were conducted on the material in a direction 90° 


to the rolling direction and parallel to the rolling direction as 
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illustrated in Figure 2. Due to the original plate thickness, the 
samples were approximately 85% of full size Charpy specimens. As 
a result the values obtained are comparative only within the group 
of samples tested. The results are listed in Table I and plotted 
in Figure 3, The heat-treated specimens tested went from a high 
value of 73 ft.lbs. (99 joules) at 70°C to a low value of 4 ft.- lbs. 
(5.4 joules) at -50°C, The as received specimens went from a 
hieh value of 47 ft, -1bs, (63. 7 joules), at.~0~- G to a low 01.6. 0) ft. -lbs, 
(8,1 joules) at -50°C, The as received specimens tested 90 degrees 
to the rolling direction (see Figure 2) went from a high value of 55 
ft.-lbs. (74.5 joules) at 70°C to a low of 18 ft.- lbs, (24. 4 joules) at 
-70°C, The two sets of specimens tested parallel to the rolling 
direction showed at sharp transition at about 10°C, however the 
specimens tested 90 degrees to the rolling direction did not show this 
transition and at no time did the energy values fall below 15 ft.- lbs, 
(20. 3 joules), This difference is attributed to the microstructure. 

At this point it may be appropriate to mention that no corre- 
lation is observed between the impact values and the fatigue life. 
Kenneford (43) also observed this, in both notched and unnotched 
specimens, 
2.5 Heat Treatment 

A number of test specimens were heat treated to increase the 


ferrite grain size. Since this particular steel contains niobium, the 
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grain structure can remain Stable at relatively high temperatures, 
Niobium can stabilize grain growth up to 1875 °F in low carbon 
steels, The test specimens were placed in annealing bags and held 


at 1750° F for fifteen minutes after the furnace reached temperature. 
After the soaking time, the furnace was shut off and the test speci- 


mens were slowly cooled down to 500 °C before they were removed 


from the furnace. The heat treatment did not remove the banding 
but the ferrite grain size increased by 17.5% from 0.022 mm to 
0.026 mm. The pearlite did not appear to be affected. 
2.6 Specimen Preparation 

Single-edge notched specimens were cut from the plate in the 
orientation shown in Figure 2. The nominal specimen dimensions were 
bel) 2u (139) (atom) slong by 13/5. (9s.52 19in) thick py 374 (19, Oornm) deep. 
The specimens were cut from the plate using a continuously cooled band 
saw. To assure the faces were all 9()° to each other, the specimens 
were finished by milling. No machining operations were performed on 
the faces perpendicular to the rolling direction. A keyhole ee notch 
used as a Stress riser was made by drilling 1/16" (1. 59 mm) diameter 
hole 3/32'' (2. 38 mm) in from the edge. The theoretical stress concentra- 
tion for this configuration is 3,25. By drilling a hole rather than 
machining a notch the orientation of the machining marks were circum- 
ferential which does not affect the stress concentration as drastically as 


does making scratches along the notch root. The metal bridge between the 
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hole and the specimen edge was cut from the inside to the outer 
surface with a jewelers saw. In this way: damage to the notch root 
was prevented, 

In order to make the crack trace easier to observe in) the 
travelling telescope, the sides of the specimens were polished down 
to No. 600 grit abrasive. The polishing was done 90° to direction 
of crack propagation. Some specimens were given a final polish 
using 6 micron diamond paste. Whereas different finishes on the 
notch root would affect fatigue life through effects on crack initiation 
time, it is felt the efiect of the two different finishes on the specimen 
Sides is negligible. 

2.7 Static and Cyclic Stress-Strain Data Determination 
2. 7.1 Static Tensile Tests for Stress-Strain Data Determination 

Tensile tests were carried out at two strain rates and at the 
four test temperatures used in this research. Two strain rates, 

0. 002 in/in/min and 0, 2 in/in/min were chosen to indicate whether 
or not the material was susceptible to strain ageing during straining. 
In addition tests at two higher temperatures, 180°C and 260°C, were 
performed at the higher strain rate. 

Photographs of the test specimens are shown in Figure 4 and 
schematically illustrated in Figure 2. The direction of testing with 
respect to the rolling direction was the same as for the fatigue tests. 


Testing was performed using an Instron testing machine which was 
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calibrated to insure that the crosshead movement could be accurately 
related to the strain. The four test temperatures were maintained 

by immersing the specimen in a suitable bath for the duration of 

each test. The two higher temperatures were maintained by enclosing 
the specimen in a furnace. Data resulting from these tests is listed 
in Table IIL 

2.7.2 Cyclic Stress-Strain Data Determination 

The cyclic tests were carried out in the fatigue machine using 
a modified specimen design as shown in Figure 4, The cyclic tests 
were carried out at the four test temperatures with the control accom- 
plished in the same manner as in the crack propagation tests. 

The test procedure involved cycling for a period of time and 
once the hysteresis loop configuration appeared to be stable a photo- 
graph was taken. Immediately after the picture was taken the eccentric 
was reset to increase the amplitude and the procedure repeated until 
the specimen failed. The data was evaluated in a method similar to 
that of Kettunen and Kocks (10), which involved plotting stress versus 
cumulative plastic strain, The results from these tests are listed 
in Table IV. The experimental data is given in Table V. 

3.0 Calculation of Stress Intensity Factor 
The stress intensity factor, K, used in the data analysis is taken 


from a derivation by Kies (ll). The form of the equation is: 


K = A M/Bh3/2fi/oe3 ~ «3 
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where A is 4,12 for cantilever loading, M is the bending moment, 

B and h are the specimen thickness and depth, and ox is equal to 

1- a/h where a is the total length of crack plus the notch, This 
equation was also used by Meyn (12) under similar conditions as in 
this program, The equation is derived for plane strain conditions 
therefore the plastic zone size should be small compared to the crack 
length and with the unbroken ligament, The validity of the equation is 
only compared with experimental data up to a/h being equal to 0.5. A 
later equation (13) which predicts the same K values for a given 
Specimen Size and crack length is accurate up to a/h =Q, 8, therefore 
the equation due to Kies is assumed accurate = the range of a/h values 
reached in this test program. 

4.0 Testing Machine and Testing Procedures 

Testing was carried out at 640 cpm on a constant deflection 
bending fatigue machine which was designed and built for this project, 
Figures 5 and 6 illustrate the apparatus schematically and as is, 
respectively. 

The specimen is mounted vertically in the bottom grips and an 
extension arm e attached to the top of the test specimen, The specimen 
is bent by a variable eccentric which activates the extension arm 
through a connecting rod which also acts as a loadcell. The drive 


shaft of the eccentric is housed in eccentric bearings which may be 
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adjusted to change the load range. The lower grip is designed to 
control the specimen temperature by either heating or cooling a 
copper extension, The test environmerit is controlled by a trans- 
parent chamber through which the desired gas is flushed during the 
test. Details of the various parts of the testing machine and pro- 
cedure are as follows: 
4,1. Eccentrics 

Two parallel eccentrics are used; one which moves the entire 
head forward or backward thus changing the center of rotation of 
the driving eccentric relative to the specimen, the other which 
changes the deflection amplitude. The amplitude of the driving 
eccentric is changed by a gear system which may be locked in place 
once the desired position is determined. 
4,2 Connecting Rod and Load Cell 

The connecting rod was machined from aluminum alloy 2024-T4 
to reduce inertia effects. At each end are self-centering bearings 
which are threaded into the aluminum rod, By threading the bearings 
into the rod it is possible to adjust the rod length, making alignment 
with the position of the extension arm simple. The connecting rod 
contains a reduced section on which two of four B, L. H. SR-4 strain 
gauges are mounted. Two of the strain gauges which serve for 
temperature compensation are mounted on a block of the same material 
which is attached to the connecting rod, 


The load cell was calibrated on an Instron testing machine as 
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outlined in Section 4. 7. 

A Tektronix 561B Oscilloscope combined with a type 3C66 
carrier amplifier were used to display the output from the load cell. 
The most sensitive settings, 10 Vom strain/div and full gain were used 
throughout the test program. The output from the scope was tied 
in with a Edin model 8002 Oscillograph which provided a visual 
recording of load versus time. The oscillograph was activated by 
a timer at intervals of 12. 5 minutes for 8 seconds in order to provide 
a load recording. 

4,3 Extension Arm 

The extension arm is made from a 3/4'' steel bar in which a 
series of various size holes are drilled to reduce the mass. A Linear 
Variable Differential Transformer (IVDT) is connected to the 
extension arm to record the deflection amplitude, and in conjunction 
with the output from the load cell a recording of the load versus 
deflection hysteresis loop is obtained. A typical load versus time and 
load versus deflection hysteresis loop is shown in Figure 7. 

4,4 Specimen Chamber 

The specimens were all tested under controlled environments 
at ambient pressure. The chamber was made of a Lucite round 
which was threaded and sealed into a base plate which is fixed to the 
lower grip. A glass viewing window and thermocouple holding device 


are designed into the chamber. A rubber boot was sealed to the top 
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of the chamber and to a plate fixed to the extension arm thus pro- 
viding a flexible seal. The impurities in the argon and hydrogen 


are shown below in parts per million: 


Argon Hydrogen 
Argon remainder --- 
Hydrogen 1 ppm remainder 
Oxygen 5 ppm 5 ppm 
Nitrogen 0 ppm So ppm 
Hydrocarbons 1 ppm 1 ppm 
Water 5 ppm 5 ppm 
Carbon Dioxide 1 ppm 1 ppm 


In addition, the argon and hydrogen were passed through a 
Matheson Model 450 dryer and a cylinder of anhydrous CaSO4 
desiccant to remove additional moisture. The hydrogen was also 
routed through a liquid nitrogen water trap. The gas entered the 
chamber through the top plate and a back pressure was maintained 
by passing the exit gas through an oil bubbler. 

4.5 Grips and ee rature Control 

In order to obtain the test temperatures the bottom half of the 
grip system was made from a copper rod encased by a steel round. 

A slot was milled out of the upper end of the copper rod and the 
Specimen was seated in this slot. The copper lobs were pressed onto 
“ie Specimen sides by set screws, thus allowing the specimen to 
contact the copper rod on three surfaces for heat transfer. The back 


edge of the test specimen was supported against a stainless steel 
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shim to reduce heat transfer to the testing machine. Another shim 
was placed at the front edge of the specimen against which two locking 
bolts pushed. When the specimen was in place the copper lobs were 
pressed by set screws up to the specimen sides, The top eight 
inches of the copper rod is centered inside of the thick walled 
cylinder by a series of set screws to minimize contact area and con- 
sequently heat losses are reduced. The copper bar is also insulated 
and sealed inside the steel holder, The. bottom twelve inches of 
copper rod can be submerged in a Dewar containing either liquid 
nitrogen, dry ice and acetone, or heated paraffin oil to obtain temp- 
eratures of -50, -10 or 71°C, The bottom grip holder is separated 
from the main frame by four small stainless steel shims, again to 
minimize heat transfer. 

The top half of the grips which is the extension arm is a box-like 
system which fits over the specimen and then is held firmly there by 
a shim and two lock bolts. 

The temperature was recorded by a copper-constantan thermo- 
couple which is held against the specimen near the crack plane. It 
was found that once the temperature reached equilibrium for the 
various liquids it remained constant. As the crack grew the heat 
transfer to the upper grip changed and the level of the coolant, or in 
the case of the higher temperature test, the input current to the 


heater, was decreased to maintain the correct temperature. A photo- 
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graph of the grips is shown in Figure 8. 
4.6 Crack Length Measurement 

The crack length was measured with a 20 power travelling 
telescope. A rule with 1/64'' divisions was placed in the chamber 
beside the specimen. With this system the crack length could be 
read to 0.0078 in (0.20 mm) with no difficulty and could be detected 
at smaller lengths. The readings were taken every 8000 cycles (or 
every 12.5 minutes) and more often when the crack growth rates 
were faster as in the cases where hydrogen was used as an environ- 
ment. 
4, 7 Calibration 

The strain gauges of the load cell were calibrated using an 
Instron testing machine over a load range of 51lbs to 100 lbs. Within 
the load range the load versus displacement plots were linear. A 
statistical analysis was carried out on the data. A mean value of 0.903 
centimeters displacement per ten pounds of load was read on the 
oscilloscope at a setting of 10 micro strain per division and full gain. 
The variance and standard deviation were 0.00077 and 0.02774 res- 
pectively. The slopes and regression coefficients for the tension and 
compression calibrations are 0, 090 and 0.9998 and 0. 0894 and 0.9997 
respectively. During the fatigue testing periodic checks were made to 


insure that the calibration did not change. 
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5. EXPERIMENTAL RESULTS 


The results of this series of tests indicate a definite difference 
in fatigue crack growth rates as the temperature and atmosphere are 
varied, Since the results were consistent within each set of tests a 
large number of repeats was considered unnecessary, however 3 to 
7 tests were performed for each testing condition. The expe aibeatal! 
data from the crack growth rate versus stress intensity relationship 
is given in Table VI. 

5) Grackulieont 

The typical appearance of the crack front is illustrated in Figure 
9. In order that the calculated stress intensity, K, be valid, the 
deviation of the crack front from linearity must be small. To obtain 
valid K values according toA.S.T.M. - E 399-70T, the crack 
trace on either surface must be equal to or greater than 90 % of the 
"average'’ crack length; ''average'’ being defined according to A. S, T. M, 
- E 399-70T point 7. 2.3 (19). Measurements were made on several 
specimens and all fell within these limits, with the average of the 
trace readings being approximately 95% of their ''average'’ crack 
length, The importance of this is that K is then fairly constant across 
the crack front since K is a function of the crack length. 


5.2 Crack Length Versus Number of Cycles 


Shown in Figure 10 is a plot of crack length, a, versus number of 
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cycles, IN, for the four test temperatures of the specimens in the as 
received condition, Figure ll shows the same plot for the heat treated 
Specimens and those tested in a hydrogen atmosphere. These curves 
can be broken into three approximate regions; an initial zone of slow 
crack growth rate, a zone of increasing crack growth rate, anda 
slower final zone as the load decreases, The crack growth rates 
(da/dN) were obtained by measuring the instantaneous slope at each 
actual data point on the crack length versus number of cycles curves, 

In constant deflection amplitude tests both the load and the 
stress intensity vary as the crack propagates thus one obtains a 
range of values for the crack propagation rates and the corresponding 
stressintensities from a single test, It is often practiced to approx- 
imate the curves of Figure 10 and ll by a straight line and correlate 
this single crack propagation rate with an average stress intensity, 
so only one pair of data is obtained from a single test, 
5.3 Crack Initiation 

Crack initiation times were determined by observation with a 
20 power travelling telescope. The time of crack initiation is taken 
as the number of cycles at which the crack first became visible. This 
was less than 0. 0078 ' (0.20 mm) of crack. For this reason all 
crack initiation times are relative to each other. Because of the 
method of defining crack initiation the relative initiation times are 


considered valid. True crack initiation occurred before this time, in 
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fact that there is a high probability that microcracks may be present 
before cycling begins. Because of the above fact normal statistics 
are not entirely valid. The crack initiation times are shown below 
and listed in Table VIL 


INITIATION TIMES IN CYCLES (X1000) 


Temperature °C Environment Cycles 

-50 Argon 818+ 9.8% 

=) Argon 58.72 15. 6% 
21 Argon 34.0 20 4 7 
71 Argon 50.0 +11. 0% 

- 50 Hydrogen 86.0 ave. of 3 tests 

- 50 Argon (heat treated) 86.0 ave. of 3 tests 
20 Hydrogen 40.0 ave. of 3 tests 
Ze Argon (heat treated) 40.0 ave. of 4 tests 


* 90% confidence 
A Rank test was used to determine whether the difference in the means 
of the crack initiation times as a function of temperature was signi- 
ficant. At the 95% confidence level it was found that the crack 
initiation times were significantly different for each of the tempera- 
tures except -10 0¢ andere Gs 

The trends indicate that 21°C has the minimum crack initiation 
time and the value increases as the temperature is decreased or 


increased. 
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5.4 Crack Growth Rate Versus OK 

Data on crack propagation rates as a function of stress inten- 
sity are commonly represented within a certain range of SK by the 
equation da/dN = C SK™, The stress intensity range, OK, is 
defined as K maximum minus Kminimum, To check this relationship 
and determine the constants C and m crack growth rates as calculated 
from the crack length versus number of cycle graphs were plotted on 
logarithmic paper against the stress intensity range, OK, Shown in 
Figures 12, 12A, 13 and 14 are the plots of da/dN versus AK for the 
20 and 70°C tests, the -50 and -10 °C tests, the heat treated Specimens, 
and for those tested in hydrogen respectively. The data is normalized 
with respect to the mean stress which varied from test to test, The 
manner in which this is accomplished is shown in Appendix l, 

The data points from each test specimen must be analysed 
separately because the density distribution of data points varies with 
4K from test to test, This would result in different weights of data 
points if data from all tests were combined for statistical analysis, In 
order to use these data properly each test must be individually statis- 
tically analysed to determine the relationship between da/dN and AK, 
Qnce this relationship for each test condition is determined the values 
of the two constants C and m can also be determined, The values of 


the constants C and m together with the correlation coefficients are 
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listed in Table VU The actual data points from each test are given 
in Table vl them value as a function of temperature 15 a 


maximum at -10°C and decreases as the temperature is increased 


or decreased, The values of C vary inversely to those of m. The 

m values for the heat treated specimens at -50 °C are high while those 
for the heat treated specimens at 21°C are in line with the as 

received condition, The effect of hydrogen is difficult to analyse in 
that only the initial crack growth rates can be correlated with the 
stress intensities before unstable crack growth occurs, Because only 
a few data points are available the values of C and m cannot be used 


comparatively for the tests in hydrogen, 


55 Fatigue Life 


The number of cycles to failure are listed in Table IX, The 
values given are the average number of cycles for the crack to grow 
three-quarters of the way across the net section, The material 
exhibited the longest fatigue life at -50 °C and the minimum at 21°C 
while the lives at -10 C and 71°C were approximately the same. The 
effect of hydrogen was to decrease the fatigue life by increasing the 
crack propagation rate. Heat treatment had little effect on the 
material at 21°C but increased crack growth at -50°C thus decreasing 


the total fatigue life. 
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5.6 Hysteresis Loops 


Load-deflection hysteresis loops were recorded during some 
of the tests. The area under these loops represents the mechanical 
energy dissipated in the specimen, the grips and connecting rod. To 
estimate the energy dissipated in the testing machine calibration 
tests were conducted using unnotched steel specimens. To calculate 
the energy dissipated in the notched fatigue specimens the area from 
the calibrated unnotched specimens was subtracted from the area 
recorded in the notched specimens, The calibration data indicated 
that approximately 80% of the total damping energy is dissipated in 
the fatigue testing machine and 20% of the absorbed energy is due to 
reversed plastic deformation, The data given in Appendix 2 will be 
used later to estimate the increase of temperature in the reversed 
plastic zone. 

5.7 Cyclic Stress-Strain Curves 

The results of cyclic stress-strain tests are given in Table IV. 
The stress versus the cumulative strain plot is shown in Figure 15, 
From the cyclic data it is observed that the order of the flow stress 
for a given strain is changed with respect to standard tensile test 


values, The flow stress for the 21°C and 71°C tests are now 


reversed in that the stress for a given strain is higher at 7] TG, hese 


data will be used later in explaining some of the observations. 
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6. ANALYSIS AND DISCUSSION 


In fatigue, cracks develop progressively from the submicro- 
scopic phase of cyclic slip and crack initiation, through a macro- 
scopicic rack propagation.pe riod, to.final fracture. 

There are two basic approaches to the development of a 
mathematical expression which describes crack propagation, one 
approach is empirical and the other is mechanistic which relates 
to the sittxucture-and, properties ofsthesmaterial... To date; the most 
successful overall approach has been based on continuum mechanics, 
However, to develop materials with high resistance to fatigue crack 
propagation the metallurgist and materials engineer must combine 
continuum mechanics with the results of fatigue studies in a mechan- 
istic approach to obtain an expression for crack propagation. 

In recognition of this, the experimental results of this thesis 
will be analysed first on the basis of continuum mechanics and secondly 
on the basis of fatigue crack propagation mechanisms taking into 
account metallurgical and environmental factors. 

6.1 Analysis Based on Fracture Mechanics 

Crack growth studies have been directed towards establishing 
a relationship between the crack growth rate da/dN and the stress 
intensity factor, K, (15), (16), or between the crack growth rate, the 
crack length, a, and the stress, Sg, (17)(18)(19)(20)(21).. The 


former relationship is represented by: 
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dafdNy="c(A Ky? (1) 
and the latter relationship by: 

da/dN = A( A 6,) ™l aP (2) 
where K = XG > Sta (stress intensity factor for a tensile loading 
of a center notched specimen), © = finite width correction factor, 
Arc=i(K maxe-cK min)» nai my and b are invariant integers, ais 
the crack length, and C and A are related to material properties. 

Frost and Dugdale (18), argued from experimental observations 
that mj = 3 andb=l1, Liu, (22) (23) used the saturation of absorbed 
hysteresis energy as a criterion for the onset of crack propagation 
andtarrivedtat my =) @and yb =i), 

A work rate model of crack propagation was used by Paris (16) 
to show that the crack propagation rate should be proportional to the 
stress intensity and after analysis of numerous experimental data 
arrived at m = 4, 

Crooker (5) has shown that the entire fatigue crack growth rate 
curve for most materials is a sigmoid when plotted on logarithmic 
coordinates as shown in Figure 16. The power-law relationship which 
holds in the central portion of the curve is bounded by upper and lower 
inflection points, The lower inflection point is taken as an indication 
of non-propagating fatigue cracks, This range may extend down to 
exceedingly low stress intensity values where the crack propagation 


rates are of the magnitude of atomic spacing in the crystal lattice. 
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(107! Q 1078 in/cyche); (10°> 21l0seemm/cycle)). The upper 

inflection point is caused by the onset of rapid unstable crack propa- 
gation prior to terminal fracture and thus places a critical range on 
the fatigue resistance of materials. 

Thus all fatigue growth rates which lead to fatigue failure fall 
within a range of stress intensity factors, bounded at high values by 
Ke, the fracture toughness, and at low values by a threshold value, 
Kth, The value of Kth is very sensitive to environment and may be 
shifted to lower values by an aggressive atmosphere. 

Recently Donahue et al, (24) modified equation(])to consider the 
threshold stress-intensity, Kth, below which no propagation occurs. 
The authors base their analysis on crack opening displacement which 
results injm stds =~ Theirvequation is of the form: 

da/dN = C(K% - Kth*) (3) 
In equation(], m and C are considered separate and independent 
quantities which represent the crack growth rate dependence on the 
stress-intensity and on material properties respectively. In the above 
case m is an invariant integer and the materials influence on crack 
growth rate resides exclusively in C. The usefulness of such an 
equation for predicting crack growth rates for a given material is 
limited because C must be evaluated experimentally. 


Originally the exponent m was thought to be invariant as 2 or 4. 


Crooker (5) has sampled a wide variety of data for steels which show 
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that the power value can vary between 1.5 and 10. However for most 
materials of structural significance the values range from approxi- 
matekyezs to’ 5: 

This shows that both m and C depend on material properties 
which must be determined by experiment, In the following it will 
be explored whether there is a relationship between these two values. 
For this purpose, data for m and ln C were taken from this study 
and the literature. The data is plotted in Figure 17...The m values 
range from 2 to 7 for a wide range of materials such as mild steels, 
high strengh steels, stainless steels, aluminum alloys, titanium 
alloys, copper alloys, and cobalt alloys. The data indicate a linear 
relationship between m and ln C of the form: 

mgs mses. AdiniG (4) 
The constant mo =-l, 3569 and the slope A = -0,1024 with a corre- 
lation coefficient of 0.993. Since the structure sensitive constant C 
is plotted in logarithmic form even small deviations can be signifi- 
cant. It is therefore proposed to modify equation(l)to read as follows: 
dafdN= C Se OK™ (3) 

In this form Sc is a structural factor which represents effects of 
metallurgical reactions and environmental effects, The value of mo 
in equation(4)may vary slightly for each class of materials or 
condition of materials. Equation(5)can be linearized to the following: 


In. da/dN =(1n.C.+.1n Sc) + mln OK (6) 
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Thus the structural factor Sc would cause a parallel shift in the 
da/dN versus MK curves. To check this postulate, data from the 
literature were evaluated. The experimental data is given in Table 
X. The measured In C values and the values as predicted from 
equation 4 are shown together with the error analysis. 

Anctil and Kila (27) investigated the effect of temperature on 
the fatigue crack propagation rate in a 4340 steel at -50°F and 80°F, 
The specimens tested at -50 °F were tempered in the range from 
400 °F to 800°F. Statistical evaluation of the data corresponding to 
equation 4 gave mo =-l. 240 and A =-0.102 with a correlation coeffi- 
cientiof-0F 99S?" =Deviation from in’ G*ranged ‘from’ -0. 293 to +@, 158 
which correspond to errors ranging from -0, 66% to +0. 365% The 
tests conducted at room temperature gave Mo =-l. 208 and A = -0,1025 
with a correlation coefficient of 0.999. Deviation from ln C ranged 
from -0, 058 to +0. 0381 which correspond to errors ranging from 
-0.173% to+0. 104%, If the test results from both temperatures are 
combined the values are mo =-l,. 212, A = -0.101 and the correlation 
Coeilicient 151059998 

Further results on 4340 steel were obtained from data reported 
by Miller (26). The specimens were austenitized at 1500°F for a half 
hour then oil quenched and tempered. The tempering temperatures 
were 200, 500, 1000 and 1400 F. The data was evaluated statistically 
for the different tempering temperatures and the results are: Mo = 


=], 504, A = =0; 1025 with a correlation coefficient of 0.999. [he errors 
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in In C range from -0. 847% to + 0. 817%. 

These results indicate that for the range of data investigated, 
tempering temperature and test temperature have little effect on 
the constant in equation 4, iie. Sc =1, 

For Sc =1, and a given mo and A, there is a fixed relation- 
ship between In da/dN and In AK as shown in Figure 18, This 
relationship is found by solving equation(4)for C and substituting 
this value into equation(]) which results in: 


lit mon 
A 

AK =(da/dN)/™ (7) 
From the above equation values of OK can be calculated for a con- 
stant crack growth rate and various values of m. Values of SK, 
and da/dN for constant m values can then be plotted, The relation- 
ship between In da/dN, 1nSK and m is expressed by a surface in 
the three dimensional plot shown in Figure 19, As noted from Figure 
18 a common intersection point for all lines occurs, This inter- 
section point is determined from the sets of data as follows: 

m-Mo 

da/dN =e (aie (8) 

Consider a constant crack growth rate and: 


mj] - m a 
1 fe) m2-Mo 


A 


ae ee 
x CAR ioe peta (9) 


All points intersect at a common value of AK. , therefore: 


ml - Mo_ m2 - Mo 
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For all the data plotted in Figure 17 and using the value of A as 
determined, the value of OK at which all curves meet is 17300 

psi Jin. The relationship below a SK of 17300 psi Jin no longer 
holds as it now predicts that materials which would show slower 
crack growth rates for a givenSK value now show higher crack 
growth rates for a givenOK, In this form no differentiation is made 
from material to material as all intersect at one value of OK, thus 
indicating that all materials have the same stress sensitivity and a 
common crack growth rate for this value of OK; This is addi- 
tional proof for the modification made to equation(], since this 
change will cause the appropriate shifts in the da/dN versus OK 
curves which are necessary to fit observed data. 

Test results of this research were also evaluated to determine 
the values of mo and A, For this purpose reference is made to 
Figure 20 which shows the crack propagation rate as a function of the 
stress intensity range for the four test temperatures. From the plot 
one observes that the lines for -50°C and -10°C intersect at one 
point while those for 21°C and 71°C intersect at another, The point 


where the 21°C and 71°C lines intersect is at a higher crack growth 


rate and lower SK value. Cyclic strain ageing should be negligible 


at the two lower temperatures as compared to the higher temperatures 
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thus one can postulate that the shift in the point of intersection is 
caused by strain ageing which takes place during fatigue cycling. 

The data for -50°C and -10°C was analysed giving Mo = -l. 228 

and A = -0.0986. The data for 21°C and 71°C gives mo = -l. 543 

and A = -0.1029, What has happened, since both slopes (A) are 
practically the same is that a parallel shift in che curves has occurred, 
The two higher temperature values have moved above the two lower 
temperature values. The shift in the two parallel lines is expressed 
by the difference in the mo values which corresponds to a difference 

in the In Cvalues, This gives Sc as in equation 5, a value of 1. 545. 
This means that strain ageing increases the basic crack growth rate by 
55%. As was previously mentioned, both C and m are affected by 
material properties, therefore superimposed to this is the effect of 
cyclic strain ageing on m, which is decreased, thus indicating a reduced 
sensitivity of the crack propagation rate to the magnitude of the stress 


intensity range. 
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6.2 Analysis Based on Fatigue Mechanisms 

Fatigue crack initiation and crack propagation is the result 
of localized reversed plastic deformation. In the case of poly- 
crystalline metals, localized deformation can occur in grains which 
are favorably oriented with respect to the maximum shear stress, 
at stresses which are substantially lower than the conventionally 
defined yield strength of the bulk material, Since grains which are 
located at the surface are constrained less than those in the interior, 
yielding at lower stresses occurs and this explains why fatigue cracks 
normally initiate at the surface. Because these grains are also in 
contact with the environment, fatigue becomes a surface and 
environment sensitive failure mechanism, Superimposed to these 
aspects is the effect of stress concentrations which result in localized 
plastic strain concentrations. Stress concentrations may be caused 
by elastic crystalline anisotropy, surface roughness, inclusions, and 
differences in elastic and plastic properties of multiphase materials. 
6.3 Crack Nucleation 

To precisely define what constitutes an initial fatigue crack is 
very difficult. Microscopic observations of the surface of specimens 
undergoing cyclic loading indicate that slip bands and microcracks 
form very early in the cyclic life, however the scale of this defor- 


mation is below the range that can be detected by conventional inspec- 


tion equipment. Therefore the definition for an ''engineering size 
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crack'"' is appropriately adopted. Manson and Hirschberg (27) 
adopted a crack length of 0. 003 in (0.076 mm) which can be detected 
with a low power microscope as the criterion for crack initiation. 

In this thesis an average crack length of 0.0078 in(0.20 mm) was 
adopted as a crack initiation criterion. 

As well, several dislocation models have been developed to 
describe the nucleation of fatigue cracks. Almost all of these models 
rely on the onset and localization of plastic deformation by slip 
leading to notching of the metal surface. Therefore, the discussion 
will be based on this ''descriptive'' model rather than the dislocation 
models, 

In single crystals of favorable orientation, cyclic loading will 
cause the formation of uniform fine slip bands on the surface and 
appreciable cyclic strain hardening occurs after only a few cycles. 
The matrix structure associated with this hardening consists of 
dislocation debris, After less than 1% o:: the total life slip begins to 
concentrate in ''persistent slip bands" in which cracks eventually 
nucleate. Persistent slip band formation is commonly observed in 
polycrystalline materials as well, The dislocation structure in these 
bands is a cell structure which is characteristic of a high amplitude 
fatigue-matrix structure. Flow stress, hardness and hysteresis loop 
measurements indicate that the persistent slip bands are ''softer" than 


the surrounding material, which means they have both a lower flow 
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stress and a lower strain hardening rate. The strain amplitude, 
temperature, and slip character influence the slip distribution and 
crack nucleation characteristics, As an example, ii the.stress 
amplitude is increased the entire specimen surface can be covered 
with persistent slip bands with a tendency in polycrystalline material 
for intergranular rather than transgranular crack initiation. The 
behaviour of persistent slip bands in commercial alloys undergoing 
low-amplitude fatigue must be similar to the bulk material undergoing 
large cyclic strains, therefore, the cyclic stress-strains response 
of alloys undergoing large cyclic strains can be assumed to be 
representative of what occurs locally at smaller strains and longer 
lives. 

The effect of grain size on fatigue crack initiation has been 
studied by several authors (28). Generally an increase in grain size 
results in a reduction of crack initiation life. When slip bands pro- 
duced by high amplitudes impinge on a grain boundary it causes an 
offset strain in the boundary which results in a high stress across the 
boundary. As cycling continues the offset and the stress increase 
until a crack initiates. The effect of grain size can be illustrated by 
assuming the same cumulative strain in two grains of different size. 
For this case the plastic strain increment per cycle O€p = A L/L 


is the same in each grain, however on this scale 2 is proportional to 


the grain size, hence the displacement (offset) Oand the resulting 
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stress will be larger in the larger grain which causes shorter crack 
initiation periods. Large differences in grain size are required to 
have a significant effect on fatigue crack initiation times, The 
increase in grain size from the heat treatment given the test specimens 
in this thesis was 17. 5% and as a result no significant difference in 
crack initiation times were observed. 

In steels strain ageing may also affect fatigue crack initiation. 
Strain ageing may be defined as the time dependent property changes 
which occur following plastic deformation, Ageing effects in steels 
result from the interaction of point defects, i.e. interstitial atoms 
such as nitrogen or carbon with dislocations. The rate of change of 
any property depends on the amount and method of deformation, 
temperature, time and the solute atom concentration, The interaction 
effects inhibit plastic flow (29). In the mortal life range, mild steels 
which exhibit distinct yield points undergo a ''softening"' at low strain 
amplitudes and a hardening at higher amplitudes, At low strain 
amplitudes softening occurs first then hardening occurs (30). The 
dislocation locking is broken down during cycling and the pronounced 
yield point disappears, Strain ageing during fatigue (dynamic strain 
ageing) depends on many variables and is rather complex. Dynamic 
strain ageing is both a cycle and time dependent phenomenon (31) (32). 
In addition it is also very sensitive to temperature, Calculations of 


the diffusion coefficient and the mean time of stay for an atom in an 
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interstitial position are given in Table XI As the temperature 
increases from 21 to 71°C a factor of 100 difference occurs in both 
the above properties. One point that is commonly neglected is the 
increase in temperature that occurs in the slip bands. If one 
considers that almost 95% of the work of deformation is converted 

to heat substantial temperature increases can occur in the slip bands, 
This may significantly modify the ageing kinetics above and beyond 
overall temperature increase effects, Measurements performed in 
this research were used to estimate the temperature increase, The 
maximum temperatures in the reversed plastic zone were calculated 
by two methods and both gave approximately the same temperature 
values, A description of the methods and the calculations are shown 
in Appendices 2 and 3, Temperature increases of at least 100 °C 
occur, and it is considered that these values are conservative in 
nature, The increased temperature effects can also be very important 
when analysing environmental effects as will be discussed later. 

The observed effects of environments on fatigue crack initia- 
tion are contradictory. Broom and Nicholson concluded that crack 
initiation was affected by environment while Brashaw and Wheeler 
concluded that the major effect was on the rate of crack propagation (33). 

In summary previous studies have shown that fatigue crack 
nucleation depends on such properties as flow stress and strain hard- 


ening and is influenced by grain size, strain ageing and possibly 
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environment. 

The results of this research will now be analysed. 

Inthe test temperature range used in this investigation, a 
wide variety of material properties occurred. Temperature depen- 
dent ageing reactions are superimposed on the effects which tempera- 
ture has on the material properties, therefore it is necessary to 
discuss these aspects before analysing other possible factors, 

The sensitivity of this steel to strain ageing has been explored. 
Figure 21 shows an increase in the yield strength and hardness during 
ageing at room temperature following plastic deformation, The upper 
yield strength increased by 23% in 100 hours while the Rockwell D 
hardness increased by 23% in 20 hours. These ageing times compare 
with an average time of 1 hour for crack initiation and 6 hours for 
total fatigue life. It should be kept in mind that localized heating can 
substantially increase the rate of ageing. At higher temperatures 
strain ageing can be so rapid that it occurs during straining leading 
to the well known blue brittleness, 

To study the effects of temperature and ageing time tensile tests 
were performed in the temperature range -50°C to 260°C at strain rates 
of 0. 2 in/in/min and 0.002 in/in/min. Shown in Figure 22 is a plot of 
yield strength and ultimate tensile strength at the two strain rates as 
a function of temperature, Also shown in Figure 23 is a plot of 


fatigue life, crack initiation, and crack propagation time as a function 
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temperature. The ultimate tensile strength exhibits the commonly 
observed trend with increasing temperatures; a monotonous decrease 
with a minimum at approximately 100°C followed by an increase in 
the blue brittle range. Because simultaneous or prior deformation 
is necessary for the ageing reaction to occur the yield strength also 
decreases as the temperatures increase. The values of yield stress 
increase as the temperature decreases below 21°C, At the lower 
Strain rates, the stress values are lower but exhibit the same trends 
indicating the material is strain rate sensitive, Although the bulk 
yield strength and tensile strength may indicate a trend similar to 
erack initiation times with respect to temperature these properties 
may not be as Significant as the flow stress since fatigue is caused by 
reversed plastic deformation, 

Flow stresses at different strains are plotted as a function of 
temperature in Figure 24, The flow stresses show additional charac- 
teristics which are different from those of yield strength and ultimate 
tensile strength. 

At the high strain rate the flow stress reaches a minimum value 
at about 80 °C and a maximum in the blue brittle range at about 190 °C, 
indicating a lowering of these critical temperatures. Strain rate has 
a Significant effect particularly at -50°C and -10°C, At these two 
temperatures the flow stress at the low strain rate is substantially 


lower than the flow stresses at the higher strain rate. The low-rate 
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flow stresses increase, at about 20 “C after which they follow parallel 
to the high strain rate values, This type of behaviour indicates 
certain ageing reactions are occurring while straining in the G Cito 


20°C temperature region. This means that two ageing reactions take 


place, oe which appears at 180 °C and the other at about 20°C. It 

is well established that the increase in the blue brittle range is caused 

by interactions of interstitial carbon with dislocations. Nitrogen has 

a higher diffusivity (see Table XI) than carbon which may cause a 

Similar increase of the flow stress at 20°C and the low strain rate. It 
is believed that nitrogen has the greatest effect on the strain ageing pro- 
cess (38) and coupled with the much greater solubility (approx. 100 times) 
(35) this lower temperature low strain rate effect is plausible, Data from 
Pugh et al, (36) may partially substantiate this observation. Tensile tests 
on high purity iron (C 0.0034%, NO. 001%) were performed at two strain 
rates as a function of temperature. A decrease in engineering fracture 
stress and tensile strength occurred at the low strain rate below room 
temperature, Although detected in the course of this work a detailed 
analysis of this aspect is beyond the scope of this thesis. 

Referring to Figure 23 which shows fatigue crack initiation life, 
and total life as a function of temperature, a minimum in fatigue crack 
initiation life occurs at 21°C. This decrease may be due to the 
decrease in yield strength or flow stress, but the further increase with 


increasing temperature is most likely due to cyclic strain ageing. 
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This is evident from the cyclic yield stress data shown in 
Pigure 25, ) The cyclic yield stress jas a function of temperature 
indicates a trend similar to crack initiation time and total life which 
is shown in Figure?23,., Asminimuny oceurs in the cyclic yield stress 


atzC. The cyclic flow stress, however, does not exhibit a pro- 


nounced minimum, The increase in the cyclic yield strength above 


21°C could result from the increased rate of strain ageing. 


A comparison of fatigue crack initiation life with other para- 
meters such as strain hardening exponent, grain size, temperature 
and environment is given on the following page and is also shown in 
Table XIL 

From this table it appears as though the greatest effect on crack 
initiation time is the test temperature, The environment and grain 
size do not have any Significant effect on crack initiation time and no 
apparent correlation is observed between the static, cumulative cyclic, 
or non-cumulative cyclic strain hardening coefficients, It is interesting 
to observe that the strain hardening exponent, n, from the static tensile 
tests all show a sli ghtly higher value at the lower strain rate which 


indicates that ageing is occurring during the test duration. 
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Ca Fatigue Crack Propagation 


Generally fatigue crack propagation is divided into two stages, 
called Stage I and Stage I], Fatigue cracks form as slip band cracks, 
(usually crystallographic unless very high amplitudes are used), in the 
manner discussed in the previous section. Stage I propagation is 
favored under conditions of low stresses or corrosion fatigue below 
the air fatigue limit (37). During crack propagation in Stage I slip 
band cracks continue to grow along planes of maximum shear stress, 
In polycrystalline materials, Stage I involves many individual slip 
band cracks which eventually link up. At this point the Stage II 
propagation begins during which the cracks grow along planes of 
maximum tensile stress (generally non -crystallographic) until it 
reaches a critical crack length for which the next load peak produces 
a tensile failure of the specimen, The fraction of the total fatigue 
life spent in Stage II crack propagation increases with increasing load 
amplitude. 

It is commonly accepted that strain hardening associated with 
localized plastic deformation are the dominant factors which control 
the rate of fatigue crack propagation (38), Laird (38) has proposed 
a plastic blunting model of fatigue crack propagation which is supported 
by many observations, This model is schematically illustrated in 
Figure 26. During the tensile part of the load cycle high localized 


plastic strains occur at the crack tip. These cause a blunting of the 
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crack tip and strain hardening in this region, Upon reversal the crack 
faces are forced together without complete reversal of the load. This 
causes the crack to advance by an amount which is proportional to 

the difference of the new surface formed by the extension in the tensile 
stroke and that regained in the compression part of the cycle, The 
full compression resharpens the crack and prepares it for the next 
tensile cycle, This 'plastic blunting" process results in the ripples 
or striations often observed on fatigue fracture surfaces, 

Such striations were observed in scanning electron fractographs 
and transmission electron fractographs as shown in Figure 27, It 
is generally agreed that each striation corresponds to one cycle of 
crack growth, thus localized crack growth rates can be calculated by 
measuring the distance between striations. Measurement of the 
striation spacing on several electron micrographs gave a crack pro- 
pagation rate of 107° in/cycle which agrees very well with the macro- 
scopically measured values. Figure 28 also shows that striations 
occur in isolated patches while the rest of the fracture surface con- 
tains dimples and irregular markings. 

The plastic deformation at the crack tip is localized within a 
small zone which is embedded in an elastic region, On increasing the 
load amplitude the size of the plastic zone increases as well as the crack 
propagation rate. Crack propagation models based on fracture mechanics 


attempt to relate the crack propagation rate to the plastic zone size by 
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elastic-plastic analysis, however, there is a distinct difference in 
the plastic zone size caused by unidirectional loading and that caused 
by fatigue loading, 

When a specimen containing a crack is subjected to zero-to- 
tension loading the plastic zone at the crack tip is subjected to con- 
stant strain amplitude cycling, Under these conditions the tensile 
mean stress present in the plastic zone on the first cycle gradually 
relaxes to zero during subsequent cycling, This is illustrated in 
Figure 29, Therefore after several cycles the plastic zone is subject 
to reversed stressing rather than zero-to-tension loading, Conse- 
quently a large zone forms on the first cycle and decreases on 


"reversed plastic zone'’, Rice (39) cal- 


subsequent cycling to the 
culated a four fold reduction in the size of the plastic zone for an ideal 
elastic-plastic material, but this analysis disregards changes due to 
strain hardening and strain ageing. The sequence of events is 
illustrated in Figure 30, On the first tensile stroke a large plastic 
zone is formed due to the high strain and stress concentration resulting 
from the sharp crack, If crack tip blunting is neglected, the stress 
concentration factor is effectively infinite and reversed plastic flow 
commences with the first increment of load reduction, creating a new 


plastic zone which is embedded in the original zone caused by the first 


tensile loading. 
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This initial or ''monotonic plastic zone" is larger at the 
surface because of plane stress conditions than in the interior of 
the crack front where plane strain conditions exist. An illustra- 
tion due to Liu (21) is shown in Figure 3l, The size of the plastic 
zone is also affected by strain hardening or softening, strain ageing, 
temperature, and the environment. For example, in steels which 
exhibit pronounced yielding plastic deformation initiates by forma- 
tion of Luders bands. Figure 32 shows a photograph of the surface 
plastic zone in a steel specimen from this series, One observes 
the propagation of Luders bands at the Aechoan ban boundary. 
Furthermore propagation of the plastic zone occurs in increments 
as indicated by the sequence of curved bands along the crack boundary. 
Figure 33 shows the surface plastic zones for the four test tempera- 
tures. Shown in Figure 34 are the reversed piace zones as revealed 
by etching which indicates the dense substructure within the plastic 
zone. A large difference between the 21°C argon environment and the 
21°C hydrogen environment is evident. The data for the four tempera- 
tures is plotted in Figure 35, The size and shape of the reversed plastic 


zone were measured and the data is listecl as follows and given in 


Table x Lit, 
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SIZE OF PLASTIC ZONES IN STEEL SPECIMENS 


Temper- Radius of Plastic Zone, in, (mm) 
Environment ature °C reversed rr surface rg roy i+ 
Argon -50 4, 5x10-4(0. 0114) 0.0251 (0. 636) 56 
aii} 7. 0x10 -4(0. 0178) 0.0270 (0. 685) 38 
2 2. 0x10 -3(0. 0510) 0. 0368 (0. 990) 18 
wal 1. 0x10 -3(0. 0259) 0. 0490 (1. 24 ) 49 
Hydrogen - 50 not detectable 0. 0335 (0. 85) 
725 " " 0.049 (1.24) 


Both the surface plastic zone and the reversed plastic zone 
increase with increasing temperature except in the temperature range 
where strain ageing is significant, i.e. at 71°C, and here the reversed 
plastic zone decreased in size. The surface plastic zone is not affected 
by strain-ageing which is plausible since it forms on the first tensile 
stroke where strain-ageing is not a factor. The surface zone is 
actually larger here and this can be related to the decreased yield 
strength at 71°C. The effect of strain-ageing on fatigue has been 
studied by several investigators (32) (40) (41) (42) (31) (29) (43). No 
complete agreement between the investigations has been reached and 
this can be well appreciated as dynamic strain-ageing is a complex 
problem related to many variables. However, Wilson (29), 1970, has 
concluded that strain-ageing does contribute to fatigue strengthening 


"at around room temperature" by limiting the spread of fatigue damage. 
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This result tends to correlate with the reversed plastic zone sizes, 
as the size at 71°C where strain ageing is significant, is reduced, 

A comparison was made between the plastic zone sizes, crack 
initiation life, and total fatigue life. No consistant correlation 
could be found between the surface zone size and the fatigue life. 
This was also observed by Weiss and Meyerson (44) who made 
measurements of the surface plastic zone. The inverse of reversed 
plastic zone size, total fatigue life and crack initiation time are 
plotted in a relative manner as a function of temperature in Figure 
36. It is evident from this figure that fatigue life and crack initia- 
tion life ‘correlate with the reversed plastic zone size, i,e., as the 
zone size increases fatigue life decreases, This can be further 
expanded to include an observation made by Irwin (45) that the size of 
the localized plastic zone at the tip of a crack is proportional to the 
stress intensity, to yield strength ratio squared. Tables III and IV 
give the static and cyclic tensile test data, From the static tensile 
data it can be observed that for a common stress intensity the surface 
plastic zone size should decrease as the temperature is decreased 


below 21°C, since the yield strength increases. However, above 21°C 


the yield strength decreases which would indicate a larger surface 
plastic zone. The static tensile data correlate with the Irwin relation- 
ship in that the surface plastic zone follows the correct trend with 


respect to yield strength. As previously mentioned the surface plastic 
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zones sizes do not correlate with the observed fatigue life. Therefore, 
it can be concluded that the static tensile data are not sufficient to 
explain the results, If the cyclic yield stress values are used in the 
relationship due to Irwin, the trend that the largest reversed plastic 
zone size occurs at 21°C for a given stress intensity is observed, It 
follows from this that as the stress intensity increases the crack 
growth rate increases as does the reversed plastic zone size. Another 
interesting point is indicated in Figure 15 which shows the static and 
cyclic stress-strain curves, The order of the curves with respect to 


temperature is modified in that the flow curve for 71 °C is above that at 
21°C for the cyclic data. This falls in line with the observed trend in 


fatigue life which can be correlated with the amount of plasticity. This 
means that for a given stress the amount of strain at 71°C is less than 
that at. 2)0G, 

The results of the above discussion definitely indicate that the 
use of cyclic data and reversed plastic zone sizes are much better suited 
to predicting observed fatigue behaviour than static tensile data and 
surface plastic zone sizes, 

It is interesting to note that hydrogen increases the surface plastic 
zone size but decreases the reversed plastic zone sizes to values which 
are below detectability. This indicates the strong interaction of 
hydrogen with the reversed plastic strain which results in the brittle 


nature of crack propagation in hydrogen, This is also visible from the 
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scanning electron micrographs shown ir Figure 37, In hydrogen 
fatigue striations are absent which indicates a low level of plastic 
ceformation, There have been at least 3000 papers (48) concerned 
with the effect of hydrogen in metals which have resulted in three 
proposed basic interactions between hydrogen and metals; however, 
most of these papers have dealt with the effect of hydrogen under 
static loading conditions or as in the case of dynamic loading,- the 
effect of hydrogen on high strength alloys. 

Elsea and Fletcher (47) have observed under static loading that 
crack propagation is a discontinuous series of crack initiations and 
propagations, Williams, (48) found no temperature dependence on 
crack growth rate between -10 and 59 °C for a 4130 steel in the presence 
of hydrogen. Bernstein (46), ina review, states that hydrogen can 
cause embrittlement between -100 and +100°% with the most severe 
effect around room temperature. Vitovec (49) has found the plastic 
deformation in the presence of hydrogen eccelerates attack, In creep 
tests the reduction in area and rupture times were also reduced, 
Jackson (50) has also found that straining in the presence of hydrogen 
can reduce the breaking stress of steel, It has also been found that 
tensile stresses across iron and steel membranes significantly increase 
the solubility of hydrogen in these membranes (51). 

Just recently Beachem (51), 1972, has proposed a new model for 


hydrogen assisted cracking. The model suggests that the presence of 
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sufficiently concentrated hydrogen dissolved in the lattice just ahead 

of the crack tip aids whatever deformation processes the micro- 
structure will allow. This new model does not support the embrittlement 
theories in that it predicts hydrogen unlocks dislocations thus allowing 
them to move or multiply at reduced stresses rather than lock them in 
place. The model also proposes that the fracture modes are dependent 
on the chemistry of the steel, the heat treatment, the microstructure, 
the crack tip stress intensity factor, and the rate of supply of hydrogen 
to the crack tip which determines the concentration there, 

With this background the test results will now be discussed, The 
crack length, a, versus number of cycles graphs are shown in Figure 
ll. The OK versus da/dN curves for the specimens tested in hydrogen 
are Shown in Figure 14. As can be observed the fatigue lives are 
considerably shorter at both -50°C and 21°C in the presence of hydrogen. 
The crack growth rate increased by an order of magnitude under these 
conditions, 

The test results indicate that an incubation period is required, 
either for the hydrogen to attain some critical concentration within the 
material or for the stress to reach some critical value (46) (52). The 
result of a 21°C test which was begun in argon then changed to hydrogen 
indicate that the presence of a sharp crack, i.e. higher stress or 
strain, reduce the incubation time required, as unstable crack propa- 


gation began very soon after the change. A photograph of this particular 
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test is shown in Figure 38. The surface plastic zone is almost non- 
existant at the start where argon was used but becomes evident when 
the hydrogen is introduced and rapid crack growth begins, A critical 


stress ntensity factor range (A Kic) of 20 - 22 ksiJin is necessary 


for the onset of rapid propagation. This stress intensity range under 
the the other test conditions did not cause rapid crack growth, 
Reference will now be made to the Table of reversed plastic zone 
sizes and surface plastic zone sizes, These are also illustrated in 
Figures 33 and 34, These figures and tables indicated that hydrogen 
increases the surface zone size but substantially decrease the reversed 
plastic zone size and the amount of fracture surface deformation, The 
surface zone size increase tends to agree with Beachem's new model 
as this zone is formed on the first tensile stroke which simulates the 
semi-static wedge loaded tests used in his report. This increase in 
surface plastic deformation may result from hydrogen reducing surface 
layers which inhibit dislocation motion, The decrease in plasticity 
as observed on the fracture surfaces and from reversed plastic zone 
measurements do not agree with Beachem's model. Under dynamic 
loading there are several conditions which would make a bad situation 
worse. These are: the presence of hydrogen while plastic deformation 
is occurring, the formation of a large plastic surface zone which pro- 


vides a greater surface area for hydrogen to diffuse through, anda 
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very important consideration which is the local temperature rise due 
to cycling which can accelerate the diffusion of hydrogen into the 


matrix ahead of the crack, 


6.5 Analysis of the Tomkins Model 


The fundamental problem in developing crack propagation theories 
is the determination of the displacements in the vicinity of the crack 
tip and relating these to the nominal stress or strain by a realistic 
stress-strain law. This problem has been approached in several 
different ways. These include the Bilby-Cottrell-Swinden (53) theory 
of continuously distributed dislocations, elastic-plastic analysis, and 
elastic-plastic measurements. Many of these theories give basic 
crack growth laws of the form da/dN = C OK™ which is identical to 
equation(l). In these laws m is a constant, not borne out by experiment 
and C is the only parameter which depends on material properties. 
From this point of view Tomkins (54) proposed a model which comes 
closest to reality. He used the Bilby-Swinden analysis, a power 
hardening law, and a plastic decohesion crack advance criterion (55). 
The power value m is a function of the cyclic strain hardening behaviour. 
For these reasons the Tomkins model will be used to analyse the data 
of this research, It is therefore expedient to review Tomkins analysis 


before using his equation to evaluate the experimental data. 
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The model is derived for a ductile metal under tension-com- 
pression loading. A schematic diagram of the crack tip configuration 
is shown in Figure 39, It was noted in some photographs (Figure 40), 
that slip bands of this nature were present, 

During the tensile half of the cycle plastic flow occurs on two 
narrow bands at + 45° to the crack tip. On increasing the strain to 
the tensile strain limit a new surface is formed by shear decohesion 
along the inner edges of the flow bands where the shear eee iene 
is amaximum, The theory assumes that on the closing half of the 
cycle,the shear flow is reversed without significant re-cohesion, It 
appears as though this assumption is justified (33). 

The Dugdale (55) model of plastic cohesive forces is applied to 
the fatigue problem since the shear strain is concentrated in a narrow 
band. The plastic zones are replaced by distributions of plastic shear 
stress, S, acting on the + 45° lines, of length D. For the plane strain 
case the exact form of S is not known, but it is assumed that it should 
be above the maximum nominal shear stress of material near the tip 
and approximately at the flow stress at the other end of the zone of 
length D, The equilibrium stress condition is given by Bilby and Swinden 
(53), from an analysis using linear dislocation arrays to represent a 
crack, The crack is represented by a continuous distribution of edge 


dislocations with their Burgers vectors perpendicular to the crack plane, 


: i ; a ; os 

pee wey 
eee i te i 3h Bia 
coinweugive> gi Apps $aF 40 mekaelb attaaoros A pee 
fle orsityith Withe sarceitg amoe wi beaten ait eae ee niyo 
JrseetG sasy eee edd vo wbaad ile dit ‘= 


Ras fo #40990 Woll vitesig Wave 2A) to ten seas oft gatsad 7 


_ , 
is sozt */ 


e 


of tietie od? qriessioat pO .qit Aosta edt oF Oap + bf ann ~ 
2 cs 
pisasoaa seeds vd bamrot 6t sositve won 6 jirotl nisada athenet oi 


ingthaty nieve: tsadé add evede ehasd wot oft to eeghs S80 ‘edt gota 


> 


aid Jo Med gatnolo odd no dedi asians yrosn ofT ccmnnikers 8 ef : 
’ 


fl moresdao-0% Ineritingte sodsiw baarevs7 vl won se£918 eid), olayo ~ | 


. A, ° 
Ee) ballitent st noliqnwesa eidt dguodt as ereeqqs _ 
bade ! 


‘ 


ov Dsilggs af exdx0i avinedas oiteale Yo (shot (84) oteingut sat 


wottan & di Baldxineonos #\ aisite tasde ol) pone meldove oegiist odd. 
i 


7 a 
= 
t rita. ie) 
aéede oiteelq 20 anditodixtaih yo beorliqet ets sence 2heslq oAT a ary 


hissieonsiq sdiuot Cl digest ito ,aonil Oep T ad? rio gatos @ site 


a 
T) 


hinorde ti reid bounweer an tf i th om tae al re raat siuencal 


: 
° i > 


qitiedt son leit sera - aga; 165d meee crn asa yee 
it Riite nite ni ts bal cite worl eas ts etosarmi 
ie etait ie aoitibes edad, cruise tliop s st a 
a ase mere rary? sisylene ms (a ow 


i 


t Hy 


. suman 8 ved battiaBorqer af wlosra ¢ i? 


bl 


The dislocations in the regions of plastic relaxation are assumed to 
experience a shear resistance acting along the slip lines and to be 
edge dislocations with their Burgers vector making an angle + © with 


the crack, The equilibrium stress condition is given by: 


Dia = JE [sec(ws /282) - 1] (12) 


Equation(12) gives the ratio of plastic shear band length, D, to crack 
length, a, for a given applied stress,6, and plastic shear stress. 
Experimental results indicate that S can be approximated by S) an 
average stress acting over the length, D, and equation(12) becomes: 

D/a = sec (m6 /282) - 1 (13) 
Equation(13)is of the form: 

A=secx-1 (14) 

If the G/ 2S ratio is small then equation(14)can be expanded by the power 
series for sec x and the result is: 

Dia= 2/8 (6/25)? [14 2/8 (6 /25)7] 05) 
Equation(15)is approximated once more in that the second term is 
multiplied by 1/32 whereas the value from the secant expansion is 5/192, 
Due to the nature of the secant expansion only the first two terms are 
significant. Up to this point the deformation zone size D, as a function 
of the crack length and stress is approximately known, A relationship 
between the amount of de-cohesion ff , and the deformation zone size is 
necessary. At the crack tip where the applied strain exceeds the fracture 


strain de-cohesion begins, implying that increased plastic strain is 
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accomodated by the formation of a newcrack surface, The amount of 
new crack surface is given by: 

Oe (16) 
where ¢ is the crack growth inerement perscycle. The fis a geometric 

factor which is approximately 1. 5 but due to the effect of the free surface 
reducing the transverse strain in the crack tip region is assumed to be 
close to |) theretore: 

See (17) 
Equation(17)is found to be generally true from experimental results, 
This is partially substantiated by Figure 36 which shows the size of 
the reversed plastic zone and fatigue life plotted versus the temperature. 
As the zone size increases the fatigue life decreases which is what 
equation(17) predicts. 

The cyclic stress-strain curve of a material can be represented 
by the following relationship: 
A 6 = B( Ret: (18) 
Since crack growth occurs in the tensile half of the cycle G = 46/2, 
this is substituted into equation(15) which gives: 

D/a =117 2/32 (A6/2S) [: +p °/32 (06 /28)2) (19) 

The assumption that the crack extension per cycle,is Se pd is now 
used, As the crack grows, the crack length, a, increases and by com- 
bining equation(17)and(19) the crack growth rate per cycle is: 


da/dN = AED =1E_(6)° a/32(28)2 C +17 2/32 (A 6/23)°) (20) 
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The value of 2S is replaced by an equivalent tensile stress T. substi- 
tution of equation(18) into equation(20) gives: 
da/dN = 12/8 (B/2T)2 0 +9°/8 (B/2T)" OE ant) ena eaencan 
Unless very high strain fatigue is considered A6/2T C€ 1 and equation 
(21) be ee 

da/dN = 1°/8 (B/2T)2 nea oe (22) 
Equation(22)can also be expressed in terms of the stress range as well 
by substituting A€,=(4 6/B)i/nt, This results in: 

da/aN = W2A62" +1/0), ajg BY (or)? (23) 
A comparison of equation(23)and equation{l)using test data of this 


investigation will be made in the following section. 


For comparison equation(23)may be rewritten to compare with 


equation(): 
da/dN=ClA6é™ta (24) 
whe re mt = (2n' +1) /n! and (24a) 
ce 2/8 BM™ (at)? (24b) 


For comparison equation) reads: 
da/dN = C (QA K)™ 
where C and m are experimentally determined. The stress intensity 
for notched bending tests used in this research is given by: 
K =A6 [i/a-a/n)? : (1-a/n)>| 1/2 (25) 
where h is the depth of the specimen. For the purpose of comparing 


results a crack length of a = 0, 2 inches and a stress of 136/70 sp81.15 
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assumed, This gives a stress intensity range of 20,000 psivin. To 


evaluate the Tomkins equation the cyclic strain hardening exponent n} 


i as calculated 


as calculated from the total accumulated strain, and n 
from the strain amplitude of each cycle were used, The data is listed 


below and is also given in Table XIV, 


Crack Propagation Rate (in / cycle) 


Test Temperature G Experimental Tomkins 
ae ty ae wearers 5 
-50 ToS S107! 7.02 x 107° 2.65 x 1077 
-10 1,90 x 10°, 1el45e107 1.90x107- 
21 eel Om 4,36x10_¢ 1.41x10_ 


is 1eSO rex lO 5, Oz x LO lesen 100 


The Tomkins equation gives crack propagation rates which are 
reasonably close to the experimentally determined ones, provided that 
cyclic stress-strain data based on accumulated strain are used, The 
trend observed with temperature is not well represented by Tomkins 
analysis but this is not surprising when considering the approxima- 
tions made in the derivation. Furthermore, strain cycling experiments 
may not fully represent the cyclic strain ageing effects of a small 
plastic zone at the tip of a crack, Ina further comparison the exponent 
m from equation(l)and the exponent m,; from equation(24a) were plotted 
as a function of temperature. These are shown in Figure 4l, The 
curves show that m; calculated from the cumulative strain hardening 


exponent n! is somewhat larger than the experimental value but the trend 
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with temperature compares favorably. Strain hardening exponents 
calculated from the strain amplitude show no correlation, 

A further check of the Tomkins equation was made by checking 
the relation between In C! and m, according to equation(4). These are 
plotted in Figure 42. One observes that the values calculated from 
the Tomkins equation fall slightly above the experimental values but 


the trend is the same and linear relationship holds, 
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SUMMARY AND CONCLUSIONS 


The main results of this investigation which pertains to cyclic 


bending fatigue tests in environments of argon and hydrogen in the 


temperature range from -50 to 71°C are summarized as follows: 


Is 


Fatigue crack initiation life and total fatigue life in argon 
decrease with increasing temperature until a minimum is 
reached at room temperature from which they increase as 
strain ageing becomes significant at higher temperatures. 
Crack initiation times are not affected by hydrogen. 

The monotonic EGE zone at the surface is very much larger 
than the reversed plastic zone. 

Hydrogen increases the size of the monotonic surface plastic 
zone, but decreases the reversed plastic zone to values which 
are below the range of microscopic observation, This concurs 
with the brittle crack propagation mode observed by scanning 
electron microscopy. 

Within a certain range of stress intensities, the crack propaga- 
tion rate can be expressed by the equation, da/dN = C(4 K)™, 
The power value m represents the sensitivity of the crack pro- 
pagation rate to the stress intensity range. 

A decrease in the test temperature causes a decrease of the 
crack propagation rate and a decrease of the sensitivity to the 


stress intensity range. 
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Cyclic strain ageing causes a decrease in the crack propaga- 


tion rate at 71°C compared with that of 219°C and a decrease in 


the sensitivity of the crack growth rate to the stress intensity 
range, 

There is a relationship between In C and m of the form mo - m = 
A 1n C which is the same for a material in different heat treated 
conditions and at different test temperatures provided no strain 
ageing occurs during cycling, 

Tomkins model for fatigue crack propagation can reasonably 
well represent experimental results but refinements of the 
mathematical analysis are possible and necessary. 

A more detailed study into the effects of localized heating in the 
crack tip area is required so that one can better understand the 
tip conditions for crack growth analysis. 

Crack eee life, total life, and crack propagation life are 
proportional to the inverse of the reversed plastic zone size. 
Mean stress has a Significant effect on fatigue crack propagation 


rates within a certain stress range. 
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TABLE I 


CHEMICAL COMPOSITION OF THE TEST MATERIAL 


ELEMENT WEIGHT % 
Carbon 0,2 
Manganese 0. 82 
Niobium O4012 
Phosphorous 0. 006 
Sulphur 0. 021 
Tin JO Oar 
Nickel 0705 
Chromium 0, (0:5= 
Molybdenum 0. 05% 
Copper 0705* 
Vanadium 0. 005% 


GIVEN MATERIAL PROPERTIES 20 °C 


PROPERTY VALUE Mega N/m¢ 
Yield Stress 48,000 psi S51,.2 
Tensile Strength 65, 000 psi | 448 4 
Fracture Strain ;205 in/in at .0L0 in/min 


Reduction of Area 33% 
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TABU Ee 


CHARPY IMPACT TESTS 


All values given are the average of two tests. 


Set #1 - Condition of material: 954°C full anneal 
Orientation: Parallel to the rolling direction. 


Temperature Energy 
sed @ Foot-pounds Joules 
70 T3550 eke) 
ae 55.0 74,5 
0 10. 0 He Foy) 
- 10 8, 8 11.90 
-25 ey 8.8 
- 50 4.0 5. 42 


Set #2 - Condition of material: as received 
Orientation: Perpendicular to the rolling direction, 


Tempe rature Energy 

C Foot-pounds Joules 

70 55 14.5 

ie 35 47.5 

0 a 39.9 

-25 20 Ziad 

- 50 fae} Slee 

- 70 18 24, 4 


Set #3 - Condition of material: as received 
Orientation: Parallel to rolling direction. 


Temperature Energy 
AG Foot-pounds Joules 
70 47, 0 G35.0 
oe 30,0 40. 7 
0 ole ee 
-25 8.-0 10,35 
- 50 6. 0 See 
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woltossth ggilier sifot valeothmaqrsT steifethaiy® ~~. 
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ture °C 


- 50 
- 10 
Zul 
Seto 


oil 


Uhr Alois Bs) 
O50 
0. 1867 
0. 2126 


Ab Lev 


Fracture Strain 


0.01862 
OF04 725 
0.01676 
0.01754 
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4, 55x10 (314) 
3, 28x10* (226) 
2) cipal 
2, 75x10~ (189) 


CYCLIC STRESS-STRAIN DATA AS A FUNCTION 
OF TEMPERATURE AT 640 C.P. M, 


4. 99x107 (344) 
4, 42x107 (305) 
4, 40x10* (304) 
4, 30x10* (297) 
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Yield Strength Ultimate Strength 
in/in mm/mm psi(mega N/ nr“) psi (mega N/m) 


TABLE V 


Experimental Data From the Cyclic Stress-Strain Tests 


Temperature 
WG 


- 50 


- 10 


Stress 


psi (mega N/m*) 


6649 ( 45, 84) 
13305 (91. 73) 
21445 (147. 86) 
27356 (188. 61) 
32543 (224, 31) 
38451 (265, 11) 
45477 (313. 55) 
46217 (318. 65) 
47141 (325. 03) 
48066 (331. 40) 
48066 (331. 40) 
499 13 (344, 14) 
46212 (318. 62) 
48988 (337..76) 
48063 (331. 38) 


7079 ( 48. 81) 
13409 ( 92, 45) 
20869 (143. 88) 
28700 (197. 88) 
32800 (226. 15) 
36531 (251. 87) 
39 134 (269. 82) 
40624 (280. 09) 
41376 (285, 28) 
40624 (280. 09) 
41003 (282. 71) 


Strain 


in/in 


0. 00150 
OR0032 1 
0. 00478 
0. 00629 
0, 00738 
0. 0085 5 
OR01093 
0.01230 
0. 01332 
0. 01469 
OL OLTsi2 
0.01657 
0.01691 
O, 01802 
0. 01862 


0.00155 
0. 00303 
0. 00477 
0. 00645 
0. 00740 
0. 00887 
0. 01062 
On0d: 156 
OF-0123 7 
Dye 
OO 13-5 


Failed prematurely due to accidental notch, 


N 
cycles @ 640 cpm 


1280 
1280 
1280 
1280 
1280 
1280 
480 
427 
427 
374 
320 
320 
214 
214 
214 
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Temperature Stress Strain N 
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NOTE: Metallographic inspection of the specimens indicated that 
many micro-cracks develop during the test. This would mean 
that the true stress would be higher than that calculated had the 
true cross sectional area been known, This is particularly true 
near the latter part of the test where high strains are encountered. 
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TABLE VI 


EXPERIMENTAL DATA FROM THE CRACK GROWTH RATE 
AND STRESS INTENSITY CALCUIATIONS 
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INITIATION TIMES IN CYCLES (X1000) 


Temperature °C Environment Cycles 

- 50 Argon Ole min ©7200" 

- 10 Argon bop tea 15.°6% 
Bas Argon 34, Oc 4% 
Te: Argon SO #0 Gr nO 

-50 Hydrogen 86..0° ave. of 3 tests 

- 10 Argon (heat treated) 86, Oe Save, ofS tests 
Zi Hydrogen 207 OR ave, O18 tests 
ZL Argon (heat treated) 40.0 ave, of 4 tests 


* 90% confidence 


TABIE, Vill 


CONSTANTS, In C and m AS DETERMINED FROM EXPERIMENTAL 


DATA 
Temperature °C a nye ‘ 
- 50 4,-49°+ 0,25 -58. B72 Cone + 2. 47434 0.992 
- 10 Ge Ot F9O60'7 -80, 61415 + 6, (OM3Z 0.961 
21 Del? Bi Oecd -63, 86812 + 4, 26183 0.981 
al 2. 66 + 0. 499 -39, 64566 + 4. 84359 0.994 


AD ergo 


PATIGUE LIFE INSCYCLES FOR THE VARIOUS TEST CONDITIONS 


Environment Temperature °C Number of Cycles (X1000) to 
failure 

Argon -50 299 -F SO 

Argon PLO Ro taleae RA Se 

Argon Zul 132) hoe 

Argon ee Zoot Lie 

Argon (heat treated) -50 240 average of 3 tests 
Argon (heat treated) raat 142 average of 3 tests 
Hydrogen -50 150 average of 3 tests 
Hydrogen on 90 average of 4 tests 


* 90% confidence level 
minimum of 4 tests 
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TABEN Ss 


Expe rimental Data Illustrating the In C and m Relationship 
With Error Analysis Included 


Anctil and Kula (25) 4340 (austinize and oil quench) Tests at-45°C 


Tempering m In C lnc Gin Gy Eeiaror 
Temperature °C actual predicted_ % 
204 3.31] -44. 27079 -44, 56357 Br, AS) CATs) -. 66 
260 2.95 -41. 10966 -41,. 03761 =. 01205 Fae ae) 
S16 2.46 -36. 14817 ~36. 23839 =, O90Z2Z -. 249 
Wed 3205") -A2. 16087 -42, 01704 . 14383 Beye 
427 3.15 -43, 15419 -42. 99647 Pages a ba io 365 


Mo =-l. 23994 
A See Oe LO 
correlation coefficient = .99818 


eG.) ae eS ae EG. SS ee ee ee eee 
ests atezco: CG 


204 2.97 4-41. 20369 Av bey Ad totem) -0, 01188 -0. 0288 
260 2.49 -36. 51809 -36. 47997 “0, 036 12 0. 104 
Su Deo = 35, 0405} -33.9 1486 0. 03465 O07 102 
37 1 Owe 9-34, 053355 -34, 11217 0.05882 -0. 173 

Mig =7lc0tel 

A ==-0, 10136 

correlation coefficient = .99991 


Miller (26) 4340 815 ( austinize and oil quench) Tests at 23 eG 


93 6.73 =78.85654 -78. 94421 =O.08767. =O. 114 
260 3,29 -45.01169 -45, 39306 -O,S813( ~=-07647, 
538 Seen = AA ec bee 24422200 0. 04856 0. 109 
760 Bete Doe Ce - 51. 83020 Onsen oe Ogee 
Wig oS OS 
A SOL O2 55 
correlation coefficient = .99979 
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Experimental Data used in Tomkins Equation (54) 


Demperature mt In c} inc din Cc! ir ror 
of Test. °C actual predicted % 
- 50 5. 85 -70. 08638 -69, 85155 0. 23483 0,335 
Ie 7. 69 =o9, 15345 -88.96664 0.18690 0. 209 
Al Tieebe -85. 49979 265243450 10: 06529 0. 0764 
Hit 6) 70 -78. 34404 -78, 68192 =0,33788 -0. 432 
Mo = -O, 87384 
A = -0, 09626 
cerrelation coefiicient =.9995 
Hemperature m In.G tae, ‘Aine G terror 
of test © actual predicted % 
- 50 4, 49 =O. 51206 -57. 77986 O. 71936 [eRe stoys 
- 10 Guss | -80, 61415 =O0— OLDS 0.09362 0.116 
Pad 5. 19 -63. 868 12 - 64, 60697 0. 73685 1.157 
re! (ES of) -39, 64566 -39, 75432 =O. 108 6b Oc 4 
mM, = -l. 38699 
A =) =, 410186 


correlation coefficient = .99926 
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SIZE OF PLASTIC ZONES IN STEEL SPECIMENS 


Temper- Radius of Plastic Zone, in (mm) 
Environment ature °C reversed ry surface rg te /rr 

Argon -50 4, 5x1074%0.0114) 0.0251 (0. 636) 56 

= i) 7. 0x10-%0, 0178) 0. 0270 (0. 685) 38 

21 2. 0x10- 0. 0510) 0. 0368 (0.990) 18 

Tal 1. Ox10- (0. 0259) 0. 0490 (1. 240) 49 
Hydrogen -50 not detectable 0, 0335 (0. 85) 
21 " " 0, 0490 (1s 24) 
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PARALLEL TO R.D. 


Exgure 2, 


Schematic drawing showing the test specimen orientations 
with respect to the rolling direction, 


Fatigue test specimen (parallel to R.D.). 
Charpy impact test specimen (parallel to R, D.}. 


Cyclic stress-strain specimens (parallel to R. D.). 
Static tensile specimens (parallel to R. D.). 
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Charpy impact test specimen (perpendicular to R. D.). 
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Figure 2A, Manganese sulphide inclusion content 
of test material (X200) 
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Figure 4, 


Test specimens used for cyclic, (upper) and 
static stress-strain data determinations. 
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Figure 6, 


Photograph of fatigue testing apparatus. 
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Fracture surface of afatigue test specimen 


which illustrates the crack front (X5. 4). 
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CRACK GROWTH RATE, da/dN, in./cycle 
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STRESS INTENSITY RANGE, dK, Ksivin. 


Figure 12, Crack propagation rate as a function of stress intensity 
range for tests in argon at 20 and 70 C. 
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Figure 12A, Crack propagation rate as a furretion of stress 
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CRACK GROWTH RATE, da/dN, in./cycle 
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CRACK GROWTH RATE, da/N, in./cycle 
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Figure 14. Crack propagation rate as a function of stress intensity 
range for specimens tested in hydrogen at -50 and 20 -sGe 
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Figure 16, 
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STRESS INTENSITY AGTOR RANGE, AK 


Schematic drawing of the crack growth rate versus 
the stress intensity range which illustrates the sig- 
moidal shape (5). 
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o Mild steel (59), Thesis \ 

e High strength steels (26)( 59)(60)( 61) 

vy Cobalt base alloys (59)(62) x 

o Copper base alloys (62) aS 

w Stainless steels (60)({ 62) \ 
Q Aluminum alloys (62) . 

& Titanium alloys (61) S 
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Figure 17. Slope m versus In C for a wide range of materials and 
testing conditions according to equation@). 


= 
» 
- 
a 7 : : _ 
es a: 


——— a 


— 
—_ 
a 
a 


a 


einedT (0?) Isare ptiM © 7 
(FONDS (9 OS) oloste dignoiie dgit @ 
f (SONG) ayolis onsd tadoD gy 


7 s 
~ 


vA . (S¢) syolla ousd 19qqoD a en 
(S8)(00) alesis eeoiniate a 4 
Zé 7 hisede evotls cnunicwtA és 
/ me eyolls enviassiT a : 
ae Sree a eee a [a 
0 Os~ Of- 0a— 08- 


' Dal ; int if 


Vs vy 


bai aininaiem Ye synet shiw « 10d D ab averse mo sqole +i eae 
Waaiis Sipe tr anibrease, soisihnes gnivest - is 


Figure ic. 


CRACK PROPAGATION RATE, INCHES PER CYCLE 


8 9 te) HM l2 13 


Predicted relationship between In(da/dN) and 
InAkK as calculated from the relationship between 
In Gund mo sown tajvigure 17. 
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Mean crack propagation rate aS a 
function of QK as found in this 
research for the four test temperatures, 
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Figure 25, Variation of the cyclic yield stress and the 
cyclic ultimate tensile stress as a function 


of temperature, 
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Figure 26. Schematic illustration of the plastic blunting 
model for fatigue crack growth (38). 
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Figure 28, Scanning electron micrograph for a 
fatigue test at 2] eCura argon showing the 
localized patches of fatigue striations. 
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Figure 29. Cycle -dependent stress relaxation of a specimen 
initiatially subjected to a mean tensile stress(63). 
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Figure 30, Schematic illustration of the stress at the 
crack tip showing the reversed and 
monotonic plastic zones (39). 
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Effect of the state of stress on the shape of 
the fatigue crack front, Transition from 
plane stress at the surface to plane strain 
in thescenter of the specimen, 
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Figure 32. Surface plastic zone showing the progres- 


sive formation of Luders bands as the 
crack propagates. 
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Figure s3 Photographs of the surface plastic zones at the fou 
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Figure 39, Model of fatigue crack advance due to 
shear on 45° stress planes (54), 
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Crack length = 4/64" 
(1. 78 mm) 


Crack length = 8/64"! 
(3, 56 mm) 
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(10, 23 mm) 


C in argon. 


ce plastic zones for increasing crack length 
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Figure 41, Power values m and mt, a function of test temperature, 
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linear relationship as shown previously in 
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APPENDIX 


The Effect of Mean Stress on Fatigue 
Crack Propagation Rates. 


Several researchers have explored the effect of mean stress 
on fatigue crack propagation rates. Brock and Schijve (64), McMillan 
and Pelloux (65), Roberts and Erdogan (66), and Walker (67) have 
proposed a general equation of the following form: 

daidN, =(G kK maxugs A Kye (1) 
The explicit form of equation(1)may vary by rearrangement of the 
variables, However, the above researchers have evaluated only data 
for aluminum and their final form of equation(])indicates that m] is 
approximately equal to mj; therefore, equation(1)can be rewritten as: 
da/dN = Cj(Kmax Qk)‘ (2) 

For the purpose of evaluating the data from this research the 
general form of equation(l)is expanded in the following manner with no 
previous assumptions regarding the power values, 


Aad G (Aol 2 am) Ne (3) 


where Kmax = AK/2+Km and AK/2Z is the stress intensity amplitude 


equal to Ka and Km is the mean stress intensity. 4A Kis defined as 
Kmax - Kmin, Equation(3)becomes: 


AidN SC OAK 2 (“item An) AK ou) 
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The stress ratio is defined as A= Ka/Km, Equation (5)then becomes: 
daydNe= Cust) 2 lee yay Ane Sct (6) 

The data of this research were converted for the effect of mean 
stress, 1,e, mean stress intensity factor according to the following 
equation: 

SIN aL [v2 + 1/ay2 O K]™ (7) 
By taking log (1 + 1/A) versus log MK for a constant crack growth rate the 
power value n was calculated using a regression analysis, The results 


for the tests in argon are as follows: 


Temperature “e n By 
- 50 Gas Og aL 
- 10 doe bt Ont 
21 0, 58 Gro on 
TT 6.80 Oo) te 


Comparison of equation(6) and(7) give the following identity: 


mj] = nm 


mM] + m2= m 


Thus one can write equation(8)as follows: 


daldNea Ge () ly Alek (8) 
whe re: 
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Given below are the modified values which consider the mean stress 


effect: 

Temperature In C ave. InAKave. mave. s* InAK s* In da/dN 
S600G - 58, 57731 10, 34455 4,495 0, 07740 0. 39744 
=102¢C -81, 50162 10. 05730 6. 81 0, 06614 0. 45907 

2G =65: 32777 10. 06434 5. 195 0. 1443 0. 60097 
AL ce -41. 07833 10, 49487 2. 66 0. 15401 0. 39786 


Shown in Figures 12 to 14 inclusive are the crack propagation rates versus 
the stress intensity range for the various testing conditions, The dashed 


border lines represent the standard deviation of the stress intensity 


and the solid line represents the mean, 


s*-the standard deviation 
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APPENDIX 2 


la, ‘J a i elise ‘ : : : 
Calculation of Temperatures Maximum in the Reve: sed Plastic Zone 


Calculated from a Cylindrical Heating Source Embedded in 


a Similar Matrix 


Approximate Calculation of Temperature Maximum within 
Reversed Plastic Zone 


In this analysis, the system is approximated by a cylindrical wire 
embedded in a matrix of the same material. The heat is generated 
within the wire and dissipated into the surrounding material. The heat 
flux per second is calculated from hysteresis loop measurements, The 
volume of the reversed plastic zone, assumed to be a cylinder is then 
divided into the heat flux and the heat per second per unit volume is then 
known for each temperature studied. The amount of energy converted to 
heat is assumed to be proportional to the radius of the reversed plastic 
zone, i.e, the larger the zone the more energy converted to heat, but 
the heat per unit volume may be less, It is further assumed that at some 
distance from the center of the reversed plastic zone the temperature is 
at the ambient. For comparative purposes this distance is assumed to 
be the radius of the surface plastic zone. The temperature distribution 
within the reversed plastic zone is not constant but since the zone is so 
minute it is assumed to be at a constant temperature. The temperature 
is then calculated for the center of the reversed plastic zone. The values 
obtained show only that a temperature rise occurs. The actual localized 
temperature on individual slip planes may be much higher but the problem 


of calculating this temperature becomes extremely complex, Following 
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is the temperature calculations. 
The solution for the temperature in this case is: 


T= q(Rg* - Rr*) + q Rg@ (In A Rg - In Rg) + To 


4k 2k 
T = maximum temperature in the reversed plastic zone. 
q = heat flux cal/ cm? sec 
k = thermal conductivity cal/sec cm °C =.12 


Rg = radius of large plastic zone (cm) 
Rr = radius of reversed plastic zone (cm) 
A = ratio of Rg/Rr 


To = ambient temperature 


as, Temperature el @ 
2784 cal/cm > sec - 50 
1789 cal/cm3 sec S110 
626 cal/cm3 sec 21 
1126 calf/cm3 sec wail 


The values of the Rg and Rr are given in Section 6.4, The temperature 


values are: 


Ambient Temperature Approximate Temperature in 
Ae Reversed Zone 
- 50 iec, -G 
tO igs ee 
21 98 °C 
ie 389 °C 


The values of the temperature calculated in this manner range from 
100 °C to 390 © while those calculated in Appendix 3, using a differen: 


approach range from approximately 100°C to 500 “C. In both cases a 
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temperature rise of approximately the same magnitude occurs. From 
this it can be concluded that a rise in temperature occurs although it 
is thought to be of a much more localized and intense nature in the 


slip bands found using the above solution. 
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ARPENDIX 3 


Maximum Lempe rature in the Reversed Plastic Zone 
as Calculated from a Slip Plane Model 


Introduction 

The concept of localized heating in the crack tip region has been 
studied by Freudenthal and Weiner (56), Williams (57) and Rice (58). 
Overall or homogeneous temperature rises have been observed by Oates 
(40) and many others, The maximum temperature rise will be cal- 
culated according to the method derived by Freudenthal and Weiner as 
it applies more directly to fatigue. 
Thermal Aspects of Fatigue 

The basic criterion is that a highly localized temperature is 
developed in the front of any active slip plane resulting from the con- 
version of work into heat. When slip proceeds by motions of dislocations 
under forces applied to the external surfaces of a crystal, the work of these 
forces, W, is transformed into potential energy of elastic distortion, We, 
and kinetic energy of the moving dislocations, Wk, and the energy, Wd, 
irreversibly dissipated in this motion by producing permanent change of 
shape. Wk is neglected and for relatively large permanent deformations 
and We << Wd which results in Wd = W. The dislocations are all con- 
sidered parallel so the process is two dimensional, It is further assumed 
that the conversion of mechanical energy into heat takes place on the slip 
plane over regions of significant atomic misfit, that is over widths of the 


dislocations, and also that the aggregates of dislocations required to 


ya veer 


ms . I coon man icles 
: ian bed rn | 
0) oral 


need vad noigsx git Aned> arti ni gis 


« non) b ym wal 1 


ma + 
. ie 


iaod bositad9t Ys 


‘3 
(82) o5tM bas (V2) eemslltiw (82) x90ioW bas stionbus 


a 

S ~~ 
i 
' 


ep 

paigQ yd bevisedo need svad seni oruiat sqerry enosasgomor tof 
ae ir? 

-[no od Hiw pair aT ie eogrris prarchicea adT #redto ya yon aad: 


be tonteW bas ladinebuer td yd bevineh borlt sens ay o8 ei 


atu: 


-Sugilsl o4 Wisenib ria ai at ye 11 


has 
ai 


5 


sugte't 1-2 aqeA ats T 
- pace 


si stuleteqeéal icles yidgid 6 Jedd ai wofratirs atead mC 7 
~soo silt mort enitiuwer enelg qife svios yns to In07} oct nt t | love J 
tnudaoelsib lo anchor yd shses07g gile aedW apad oni Ax0w Yo. olerey 
seer) to Arow st {stevao s lo soosliye lenterxs ect a} bsilqqn es0t0 sbay 
(OW Hotlvei sib Nestle to yqrsnie Lettinetog ola bamrolensay ei aw e. 2 8 
DW  yyprene ont bre fw ancitecolaih gnivoen adi to yarea. 


he egiteds Inonan sq gnisoborg ‘ec nego "ahs ni beiagivetb yldne rv 


ods ATT OME Ineneeri eq ofr! ytoviaater wot bas betoalgon ns W .oqade 


‘tos ils avs ameitenoleih sdT- .W < BW ai bee aise bw > > ow 
oh? 


- 


bomuaas tofdau? shal disnatansentth owd ai seenong on aa ive 
qila sd) no soslq tesa Jade adel yprans taoinadt 
As ertthtw aevo 6r near 
(obese nots ‘ we 


141 


produce the total slip are separated by the width of an individual dis- 
location so dislocations form a piled up group. The trail of dislocations 
on a given slip plane represents a moving band heat source. The heat 
released over the band source may be estimated as follows: 

Wd = W = Taba 
where TJ is the component of applied stress on the slip plane, a is the 
area swept out by the dislocation, b is the unit step produced by the 
motion of a single dislocation and n the number of slip steps to produce 
the total slip distance ra = 11D; 

The dislocations move with an average velocity, V, and the work 
is completed during the time t = a/V so the average rate at which work 
is done W/t ~ dW/dt = TdV. The width of the dislocation in a slip 
plane is equal to mb, the length A of a moving heat source is A = nib 
and dW/dt is dissipated uniformly over this length, the uniform heat-flux 
density q released over a band of length A is q = 71 §V/nmb = TV/m. 
The velocity, V, of the moving trail of dislocations is assumed to be 
between 0.1 and 0.9 the velocity of sound in the material and the distance 
m between consecutive dislocations is roughly between 3b and 7b, The 
solution for a moving heat source in an infinite medium can apply to a 
tiny reversed plastic zone at the tip of a fatigue crack. Considering the 
band heat source of length, A , which is taken as the radius of the 
reversed plastic zone and the flux energy q moving in an infinite medium 


the heat conduction is governed by the following equation: 


salle att sorwes tao based ree lq 
vewortion an badariites aud vam serve. basd orl S r9v0 ee iss. 
? ° = i 
Saale 2W = sw 
rae i ee Le) ! 
ant «ie saelg qite acd no Reotie Boniaas to swonogenes wel Fae . wv 
ot yo Doowhorg qete finn offi ai d otseooleib ‘i vid 0 4 
bes coy ee i ‘f 
" ouboxg | od ages qile to redmnue adi a bas noitaoo feib gigaia # 10 noi a : 
1s 3 7 » 
dn > $ jometie? sh rics 
i 7 - 
lrow edd bas ,V ,\ysioolev aga7eve As saatiok avom eaotteaofeib ras _ 


fp . : 


, 
a 
erow dotdwy te 9361 Sygerevs oj of is = 4 ami) a¢lit anitnb Sotalgng eo ni 


qile ai aparece oh) lo Aibiw sdee Var = = Ib\Wh = 1\Ww Brie 
n tes 7 


arc = h, #i s5%u08 ar gnivorn & lo kh dignel!l ed) wd og {aupe ai on iq 


‘| 
it \¥P = dean\VOT 8 ai djgned jo rea # 1840 0 baewalen ie | 


sd od berwuees ai anolisooleth to ex grivom at wo Vv thot aT 
: aos Ge 
suet od) Gas lair sison odd ai bavow to ydiaolev add 0.0 bas I On 


- : o - ae ’ en LD 


at dT bas dl neewiod vider Bi sigiianiptagtierink ‘eviwo9en02 A 


4 _ A» 
6 6) yiqgs ne> muibers stiniin) as ai comuon taon Bnivory' a 
: ‘4 it Ay rr 


atlt gait. obienoD bier suigits # be qid odd ja anos shtasiq foe 
Cf _ 7 pelee 


sso, suites ous es pode wi faite is, toe soruos te 


soil 1 be 
= [oe = i 


Rascliahar sito bs mm parte pear aa acid Kas 
7 - > : : 


ae: 


142 


Kk (0% T/dx* +02 T/dy2)-vadT/dx=0 

The boundary conditions are: 
Fhe OD/dy)is a/2 for 0<xa Awe 0 and 
(eT/dy)=0 for x<0, xs sy =0 and 


limotfas liga T L540 


xX—tos y—> too 
T = 0 is taken as the temperature of the undisturbed temperature of the 
medium, Heat conduction in the direction of the medium is neglected and 
inelmaximumitemperature-T occurring at x = A , y= 0 1s with the use 
of Qiao Vira, T= T/mi(k A Via ee 
The calculation of the maximum temperature is made substituting 
in various values of A and corresponding values of TY for a given strain 
for each of the four fatigue-test temperatures. The maximum temperatures 
calculated serve only to show that the tip temperature is increased and 
the values are not to be considered exact. There are many possible 
reactions and effects that increased temperature may promote. Several 


of these are listed below and some could become a research topic: 


slp thermal softening 


ies) 


local precipitation 


oy heat straining and oxide formation 

4, increased rate of strain-ageing 

D resolution of fine precipitates under action of shear and temperature 
6. decrease in local tensile and yield strength depending on actual 


temperature rise. 
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ie Possibility of over ageing 
8, Much faster diffusion rates of interstitial solute atoms 
oe In the case of hydrogen, a much faster and deeper penetration 


into the specimen due to faster diffusion rates. 


10. Formation of compounds such as nitrides, carbides, possible 
hydrides, which required additional energy in the form of heat. 


le In the case of local heat flashes on individual slip planes, very 
severe thermal stress gradients which may be in the order of 
the tensile strength which may be a plausible mechanism of crack 
initiation, 


Galeaulation of Demperature at Grack Ti 
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moving heat source q 


Solution for Moving Heat Source in an Infinite Medium 
T= T/mK (k Av/n )1!? 


T =maximum temperature °C 


T = stress =1/2 cyclic stress at 0, 01 strain 


K = thermal conductivity 0.12 cal/sec cm a 
Z. 

k = thermal diffusivity = Ko e ="0,1388: crnfa) sec 

m = 5 


@ = density 7. 86 arigiem 
Ge =speciic heat @pliscal/ gm ~ © 


= velocity of dislocation 0,1 Vo = 0.79 Vo 
Vo = speed of sound in material = 51.3 x 10° cm/sec 
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A 


* = length of band heat source 
is taken as radius of reversed plastic zone as measured for 
each temperature 


Temperature °C 4 
750 nee Oe ie wera 
a0 1778x1072 cm 
oy BOS) S105 2 vem 
wl 2. 539 x 1072 cm 
Vion ast ex 106 cm/sec 
tal, Cpe 51,3 x 107 cm/sec 
.45 Vo = 23.1x10% cm/sec 
9 Vo = 46.2 x 10 cm/sec 
Temperature 2G a psi T (dynes/cm®¢) 
50) 22, 750 1,568 x10 
210 19, 500 l. 3444x 10 
oe 18, 750 1.292 x10 
ua 179250 1 189 x 10 
Temperature °C Maximum 


Dislocation Velocity 


Tempe rature °C 
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